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1 Introduction 


One of the interesting featnres of SU(3) flavonr symmetry breaking in the light nentral 
meson sector is the mixing of singlet and octet states in the psendoscalar mesons. The 
physical r| and r|' particles can be expressed as an admixtnre of light qnark and strange 
quark flavour eigenstates (^, 


V lb') ) V 


-sin0p \ / Ipq) \ 
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where the quark states are defined as 


Ibq) 


^^|uu + dd) and |ris) = |ss), 
V2 


( 1 ) 

( 2 ) 


and 0p is the mixing angle between light and strange quark states. In principle all 
possible SU(3) flavour singlet states that contribute to the physical particles, including 
the gluonic wavefunction, \gg)., and the heavier quarkonia wavefunctions, |cc) and |bb), 
must be considered. Of these states, only \gg) is considered to contribute to the mass 
eigenstates [2|^, and this only occurs in the heavier p' meson [^|^. Introducing the gluon 
mixing angle 0 g, the p' wavefunction becomes 


Ip') = cos0Gsin0p|pq) + cos0gcos0p|Ps) +sin0G|5'5') • 


( 3 ) 


Many phenomenological and experimental studies have been carried out to determine the 

38 - 46° |6}|n 


values of the mixing angles, and measurements are in the range 


Refs. 15-17 


also contain measurements of the gluonic mixing angle, which are consistent 
with zero, albeit with large uncertainties. Most recently, the mixing angles have been 
measured by the LHCb experiment using 3°^^ —)■ J/i[)p(') decays, and are found to be 
0p = (43.5;^!)° and 0 g = (0 ± 25)° [^. 

One consequence of the mixing is the difference in branching fractions for b-hadron 
decays to hnal states containing p and p' mesons. The gluonic mixing alters the branching 
fraction for the decays to p' mesons compared with the equivalent decay to p mesons. 
Possible Feynman diagrams for amplitudes contributing to the A° —)■ Ap decays are shown 
in Fig. with Fig. 1(a) showing the dominant b—s transition via a penguin diagram. 
Due to the gluonic contribution of the wavefunction, extra Feynman diagrams of similar 
amplitude are available for the p' decay |19 . These include the possible non-spectator 


contribution, where the light quark radiates a gluon which can form an p' meson, as shown 
Fig. |l(b) and the anomalous contribution, where an excited gluon in the decay can 


m 


radiate an p' meson, as shown in Fig. |l(c)[ The amplitudes for these extra processes 
can interfere to enhance or reduce the branching fraction to p' mesons. For example, the 
branching fraction for the B°—)■ K'^p' decajj^has been measured to be (6.6 ± 0.4) x 10 
which is over 50 times larger than the branching fraction for the —)■ K°p decay . 
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Charge-conjugation is implied throughout. 
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(a) Dominant electroweak loop diagram. 


(b) Non-spectator diagram. 


AS 
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(c) Anomalous diagram. 

Figure 1: Feynman diagrams for the AS—?- Ap^'^ decay. The non-spectator and anomalous 
diagrams are available only to the r)' through the gluonic contribution. 


By measuring the relative branching fractions of many different decays to hnal states 
containing rj and rj' mesons, it is possible to extract a measurement of the mixing angle (j)p. 
Decays of b-baryons to hnal states containing r] or r{ mesons have not yet been observed; 
however, the branching fractions of the A^ decays have been estimated to be in the range 
(1.8 —19.0) X 10“® 21 , depending on the model used to calculate the hadronic form factors. 


The interference between the anomalous and non-spectator contributions cancel in such 
a way that the branching fractions for the A° —)• Ap' and A^ —)■ Ap decays are expected 
to be similar. Using the QCD sum rules approach to calculate the hadronic form factors, 
the branching fractions are predicted to be in the range (6.0 — 19.0) x 10“®. The pole 
model approach predicts smaller branching fractions, in the range (1.8 — 4.5) x 10“®. For 
comparison, if the anomalous contribution is neglected, the branching fraction for the 
A°—)■ Ap' decay would increase to (33 — 40) x 10“® 
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This paper describes the search for the A{( —)■ Ap' and A° —>■ Ap decays and measurement 


of the relative branching fractions with respect to the —)■ K%' decay, using the 3 fb 

of data in pp collisions collected in 2011 and 2012 by the LHCb experiment. 
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2 LHCb detector, trigger and event simulation 

The LHCb detector |22p^ is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < p < 5, designed for the study of particles containing b or c quarks. The detector 
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includes a high-precision tracking system consisting of a silicon-strip vertex detector 
surrounding the pp interaction region j^, a large-area silicon-strip detector located 
upstream of a dipole magnet with a bending power of about 4Tm, and three stations of 
silicon-strip detectors and straw drift tubes 25 placed downstream of the magnet. The 
tracking system provides a measurement of momentum, p, of charged particles with a 
relative uncertainty that varies from 0.5% at low momentum to 1.0% at 200GeV/c. The 
minimum distance of a track to a primary vertex, the impact parameter, is measured 
with a resolution of (15 -|- 29/pt) M-ni, where px is the component of the momentum 
transverse to the beam, in GeV/c. Different types of charged hadrons are distinguished 
using information from two ring-imaging Cherenkov (RICH) detectors 26 . Photons, 


electrons and hadrons are identihed by a calorimeter system consisting of scintillating- 
pad and preshower detectors, an electromagnetic calorimeter and a hadronic calorimeter. 
Muons are identihed by a system composed of alternating layers of iron and multiwire 
proportional chambers 


27 


The trigger 28 consists of a hardware stage, based on information from the calorimeter 
and muon systems, followed by a software stage, in which all charged particles with 
Pt > 500 (300) MeV/c are reconstructed for 2011 (2012) data. The software trigger 
requires a two-, three- or four-track secondary vertex with a signihcant displacement 
from any primary pp interaction vertex (PV). At least one charged particle must have a 
transverse momentum px > 1.7 GeV/c and be inconsistent with originating from a PV. A 
multivariate algorithm 29 is used for the identihcation of secondary vertices consistent 
with the decay of a b-hadron. 

In the simulation, pp collisions are generated using Pythia with a specihc LHGb 
conhguration 31 . Decays of hadronic particles are described by EvtGen 32 , in which 
hnal-state radiation is generated using Photos I^. The interaction of the generated 


particles with the detector, and its response, are implemented using the Geant4 toolkit 
as described in Ref. 
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3 Event selection 

Candidate signal events are identihed by reconstructing the A{( —)■ Ap' and A{( —)■ Ap decays. 
The decay B° —)• Kgp' is used as a normalisation channel for the measurement of the 
branching fractions of the signal decays. 

The long-lived Kg and A particles are reconstructed through the Kg —)■ 7f*"7r~ and 
A—)■ p7r“ decays. The reconstruction of these particles is labelled according to where in the 
LHGb detector the decay occurs. If the particle decay products produce hits in the vertex 
detector then the candidate is classihed as Long (L)] otherwise the candidate is referred 
to as Downstream (D). Since the track resolution is different for the two categories, the 
selection is optimised separately for L and D candidates. 

The candidate p' mesons in the A{(—Ap' decay are reconstructed as p'—?■ tDtC') and 
p'—)■ 7f''7r“p, where the p meson decays into two photons. The candidate p mesons used 
in the reconstruction of the A{(—)■ Ap decay are reconstructed with the p—)■ TDn~n^ decay. 
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with the 7t° meson decaying into two photons. Since the charged tracks in the r] —)■ 7i^n~nP 
decay can be used for the trigger selection, the analysis is based on this mode rather than 
the neutral decays of the rj meson, which would have a lower efficiency. 

Events that pass the trigger selection are snbject to fnrther requirements consisting 
of kinematic, particle identification and multivariate selections. A minimum requirement 
is placed on the ht qnality of the reconstrncted vertices. When reconstrncting the B° 
and A° candidates, the direction angl^ and the impact parameter are required to be 
consistent with the reconstrncted particle originating from a PV, and the calcnlated 
lifetime should be signihcantly different from zero. The (A°) candidates are required to 
have pt > 1.5 (1.0) GeV/c. The Kg mesons are required to have px > 1-2 GeV/c and a flight 
distance signihcantly different from zero. The invariant mass of the 7f''7r~ pair should be 


within 14MeV/c^ of the known K*’ mass 20 for events in the L category, and be within 
23MeV/c^ for events in the D category, where the mass window is chosen to be three times 
the resolntion of the invariant mass. The A baryon shonld have px > 1 GeV/c and an 


invariant mass within 15MeV/c^ of the known A mass 20 for the L category, and within 
20MeV/c^ for the D category. The ri' and r| mesons from B'^ and A{( decays are required 
to have px > 2 GeV/c. Reconstructed 7r+7r~7 and 7t+7r~ri candidates are required to have 
an invariant mass in the range [0.9,1.05] GeV/c^ and candidates shonld have an 


invariant mass in a 150MeV/c^ window aronnd the known q meson mass 20 . Charged 


tracks are required to be of good quality with px > 300MeV/c. The q mesons from q' 
decays, and 7X^ mesons from q decays are reconstrncted from two photons which can be 
resolved in the calorimeter, with an invariant mass in a 50MeV/c^ window aronnd the q 
and the q or 7r° mesons are reqnired to have px > 200MeV/c. Finally, 


rO 


20 


or 7t mass 

photons are required to have a transverse energy > 200 MeV, and a large p-value for 
the single photon hypothesis, in order to reject background from misidentified 7T° mesons 
where both decay photons form a single merged cluster in the calorimeter. 

To improve the resolntion of the reconstrncted invariant mass, the fnll decay chain is 
refitted, where the tracks and displaced vertices are constrained, with the position of the 
PV of the B*’ or A° candidate hxed to the PV rehtted using only tracks not associated to 
the b-hadron decays, and the invariant mass of the Kg, A, q' and q particles hxed to their 
known masses 36 . Candidates with a poor qnality in this constrained ht are rejected. 


which removes approximately 90% of backgronnd from the L category and 20% of the 
background from the D category. 

Information from the RICH detectors and the LHCb calorimeter system is nsed in 
a nenral network to constrnct a probability that a track is a pion, a kaon or a proton. 
This particle identihcation (PID) information is nsed to rednce the backgronnd from 
misidentihed kaons, protons and pions. 


A boosted decision tree (BDT) [W 38 is used to obtain further discrimination between 
signal and background events. A different BDT is trained for each signal channel and the 
normalisation channel, separately for the L and D categories and the 2011 and 2012 data 
samples. Samples of simnlated events are nsed to model the signal decays, and events 


^The direction angle is defined as the angle between the momentum vector of the particle and the 
vector between the PV and the decay vertex. 
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in data with a reconstructed invariant mass greater than lOOMeV/c^ above the or 
are used to model the combinatorial background. 


mass 


20 


In this analysis, the numbers of candidates in both the signal and the background 
samples are limited and the performance of the BDT is improved in two ways. In the 
simulated signal sample no requirements are placed on the trigger selection. Instead, the 
kinematic distributions of the candidates are reweighted in order to match the harder 
distributions of candidates which pass the trigger selection. In addition, the signal and 
background samples are split in two, and two BDTs are trained for each channel. The 
hrst half of the sample is used to train a BDT which is applied to the second half of the 
data, and vice versa for the second BDT. This way, all data are available in the training 
of the BDTs, while any bias in the multivariate selection is avoided. 

The BDT uses a set of variables with discriminating power between signal and back¬ 
ground, including kinematic variables and vertex and track quality variables, and combines 
them into one variable which separates signal and background well. The selection is 
optimised using the hgure of merit eMVA/(f + V^b) 39 , where ejvivA is the selection 


efficiency for a particular BDT selection, a = 3 is the targeted signal signihcance in stan¬ 
dard deviations, and is the number of combinatorial background events reconstructed 
within the signal region and passing the BDT selection, found by performing an unbinned 
extended maximum likelihood £t to the sidebands of the data and extrapolating this £t 
into the signal region. 

The hnal stage is to ensure that each event passing the selection contains exactly one 
B° or candidate. Due to the poor invariant mass resolution of the mesons, it is 
possible for low energy photons from the underlying event to be reconstructed with signal 
pions to create a second candidate. This happens in about 10% — 20% of events which 
pass the selection. In this situation the candidate with the highest photon is kept, and 
all other candidates in the event are rejected. 

The efficiency of this selection is (0.032±0.001)% in the L category and (0.030±0.001)% 
in the D category. This includes the efficiency of reconstructing and selecting events which 
are simulated within the LHCb detector acceptance. 


4 Results 

4.1 Fit results and signal yields 

An unbinned extended maximum likelihood £t to the candidate b-hadron mass spectrum 
is performed on the data which pass the selection. The model used for the ht to both 
the signal and normalisation channels consists of an exponential function to describe the 
combinatorial background, and a sum of two Gaussian functions with a common mean to 
describe the signal. The ratio of the resolutions of the two Gaussian functions and the 
ratio of the signal yields are obtained from a ht to the mass distribution in simulated signal 
samples. These parameter values are used in the ht to the data, and only the resolution and 
signal yield of the hrst Gaussian function and, for the normalisation channel, the common 
mean of the two Gaussian functions are allowed to hoat. The resolution of the reconstructed 
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Figure 2: Mass distribution of the selected )• Kgr|' (ri'—)• 7t'''7r“7) candidates in the 2011 
and 2012 data, reconstructed in the L (left) and D (right) categories. The results of the fit, as 
described in the text, are overlaid. 


KgTi' invariant mass is 30.0 ± 1.5MeV/c^ for L candidates and 29.4 ± 1.3MeV/c^ for D 
candidates. For the signal channels, the resolution of the reconstructed Arj' invariant mass 
is 29.1 ± 1.8MeV/c^ (31.1 ± 4.5MeV/c^) for candidates reconstrncted with rj'—)■ 7t+7t“7 in 
the L{D) category, and 47.8±9.2 MeV/c^ (56.6± 10.9 MeV/c^) for candidates reconstructed 
with 7f^7r“ri in the L [D) category. Reconstructed Ar| candidates have an invariant 
mass resolution of 49.4 ± 3.4MeV/c^ in the L category and 47.6 ± 9.2MeV/c^ in the D 
category. The parameters of the ht are found to be consistent between the 2011 and 2012 
data samples, and so the two samples are added together to perform an overall £t. 

The possible presence of physics backgrounds has been investigated. The most likely 
backgrounds are: b-hadron decay modes to mesons with open charm and an meson, 
with a 7r° meson which is not reconstructed; the nonresonant decays to Kg or A particles 
with two charged pions which are combined with a combinatorial photon, 7r° or p meson to 
form an p^'^ meson candidate; or, in the case of the p'—Tt+Ti'y decays, the nonresonant 

—)■ Kg7r+7r~y or A{( —)■ A7r+7r“y decays. These backgrounds are rejected well by the 
BDT selection, and there are expected to be fewer than one candidate from each category 
passing the selection. 

For the normalisation channel, the mass distribution of the selected Kgp' candidates is 
shown in Fig. with the result of the ht superimposed. The signal yields are 188 ± 16 T 
candidates and 149 ± 14 R candidates, and the signal to background ratio is 1.8 in the L 
category, and 1.7 in the D category. Fig. shows the invariant mass distribution of the 
reconstructed p' mesons for these decays. The distribution is htted with two Crystal Ball 
functions 40 , and the parameters are found to be consistent with hts to the simulated 
samples, with a core resolution of 16 ± 6 MeV/c^ for L candidates and 13 ± 1 MeV/c^ for D 
candidates. 

For the signal channel. Fig. shows the Ap' invariant mass distribution for candidates 
reconstrncted in the A{( —)■ Ap' (p' —7r+7r~7) decay, while Fig.|^shows the same distribution 
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Figure 3: Mass distribution of the reconstructed r|' meson for selected )• Kgr|' (r|'—)■ 7f''7T“7) 
candidates, in the L (left) and D (right) categories. The results of the fit, as described in the 
text, are overlaid. 




Figure 4: Mass distribution of the selected A{(—>■ Ar|' (r|'—)> 7f^7T“7) candidates in the 2011 and 
2012 data, reconstructed in the L (left) and D (right) categories. The results of the fit, as 
described in the text, are overlaid. 


for candidates in the A°—)■ Arj' (r|^—t 7r+7r~ri) channel, and Fig. shows the Ar] invariant 
mass distribution for A° —)■ Arj (r| —>■ 7t+7r~7r'^) candidates. An unbinned extended maximum 
likelihood ht is performed using the same model as for the decay, with an exponential 
function to describe the combinatorial background and a sum of two Gaussian functions to 
model the signal; all parameters are hxed to the values found from hts to the simulation, 
and only the numbers of signal and background events are allowed to float in the £t. No 
signihcant signal is observed above the expected background for the A{( —Ar|' channel, 
and so an upper limit is placed on the ratio of branching fractions. From the hts, the 
combined signal yields are 1.0 ± 4.4 A°—;• Ap' (r|'—?■ n^n~'y) candidates and —4.2 ± 2.3 
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Figure 5: Mass distribution of the selected A°—)• Arj' (rj'—^ 7T+7t“r|) candidates in the 2011 and 
2012 data, reconstructed in the L (left) and D (right) categories. The results of the fit, as 
described in the text, are overlaid. 




Figure 6: Mass distribution of the selected A°—)• Arj (rj—)- 7t+7T“7T°) candidates in the 2011 and 
2012 data, reconstructed in the L (left) and D (right) categories. The results of the fit, as 
described in the text, are overlaid. 


A°—)■ Arj' (r|'—?■ tt+tt rj) candidates. Evidence is seen for the presence of the A°—)■ Ap 
decay with a signal yield of 5.3 ± 3.8 candidates, and a significance of 3.0cr. 
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Table 1: Components contributing to the scale factor a. 



+ 

1 

Ap'(7t+7r p) 

Ap(7r+7r 71°) 

Factor 

L 

D 

L 

D 

L D 

ftot(B°)/ etot(Ab) 

2.3 ±0.1 1.55 ±0.08 

7.4 ±0.6 

9.5 ±1.4 

4.6 ±0.3 3.4 ±0.2 

/b/ fhi 

2.5 ±0.2 

2.5 ±0.2 

2.5 ±0.2 

l/C, 

1 (fixed) 

0.95 ±0.04 

1.13 ±0.04 

0.5.B(Kg)/B(A) 

0.541 ± 0.004 

0.541 ±0.004 

0.541 ± 0.004 

B(ii')/B(t,(')) 

1 (fixed) 

1.71 ±0.05 

1.31 ±0.03 

a 

3.1 ±0.3 

2.1 ±0.2 

17.7 ±2.3 22.8 ±4.0 

9.5 ±1.2 7.0 ±0.8 


4.2 Branching fraction measurement 

The ratio of branching fractions can be measured for each signal decay with respect to the 
normalisation channel using 


b(b“^kV) 

jVs(Ag) t«(B°) /b 1 B(ti') 0.5xB(K°) 

A^slB") e„(Ag) A; C, e(ill'>) B(A) 


= a X 


A's(Ag) 

iVs(B»)^ 


where; iVs is the number of signal events determined from the fits to data; ctot is 
the total efficiency, which is the product of the detector acceptance, reconstruction and 
selection efficiencies; /b//ao is the ratio of B° to A° production fractions, previously 
measured by LHCb [^; l/C*^ is a correction factor applied to account for the photon re¬ 
construction efficiency in the A°—)■ Ap' (p'—?■ Tr+Tfp) and A°—)■ Ap (p—)■ T&n~n^) decays, 
which have more photons in the decay than the normalisation channel; and i3(p')/i3(p^'i) 
and i3(Kg)/i3(A) are the ratios of p' to p^'i branching fractions and Kg to A branching 
fractions, respectively, from Ref. 20 , where the factor 0.5 accounts for the fact that only 


half the mesons decay as Kg mesons. Each of these factors have been measured, as 
described in Sec. 5, and are given in Table along with the calculated value of a for each 
signal channel. 


5 Systematic uncertainties 


A summary of systematic uncertainties is given in Table The largest systematic 
uncertainty for this analysis is due to the limited knowledge of the ratio of production 
fractions, /b//ao, which was measured in Ref. 41 as a function of the A{^ pseudorapidity. 
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Table 2: Fractional systematic uncertainties, in percent, on the ratio of branching fractions. 


Source 

Ari'(71+71 y) 

L (%) D (%) 

Ari'(7r+7r ri) 

L (%) D (%) 

Ari(7t+7r T\P) 

L (%) D (%) 

B (1/0) 


0.73 

0.73 


0.73 

B (t|<'') 


0.0 

2.9 


2.3 

/b//aO 


8.0 

8.0 


8.0 

Fit model 

1.7 

0.3 

1.7 0.3 

1.7 

0.3 

Ratio of Cacc 


2.0 

1.8 


1.9 

Ratio of esei 

5.9 

5.3 

8.6 14.5 

7.5 

5.9 

^trig 


1.0 

1.0 


1.0 



- 

4.2 


4.6 

FID 

2.1 

0.9 

2.0 1.0 

2.0 

1.0 

Multiple Cand 

1.4 

1.9 

2.2 1.4 

2.0 

1.7 

Total 

10.6 

10.1 

13.4 17.5 

12.7 

11.5 


For the average A° pseudorapidity in our signal sample, the ratio /b//ao has a value of 
2.5 ±0.2. 

The systematic uncertainty on the ratio of the branching fractions of the Kg, A and 
decays is calculated from the average values from Ref. 


20 


There is an uncertainty on the number of signal decays due to the ht model. To 
evaluate this uncertainty, the parameters that are hxed in the extended unbinned maximum 
likelihood ht are varied within the uncertainties obtained from the ht to the simulated 
samples, and the systematic uncertainty is the relative change in the yield obtained from 
the ht. This is found to be a small effect, 1.7% for the L model and 0.3% for the D model. 

The selection efficiencies are calculated with independent simulated samples, produced 
with diherent trigger conditions. The uncertainty on the measured ratio of efficiencies is 
the statistical uncertainty due to the number of simulation events generated. In addition, 
there is a systematic uncertainty due to the measured efficiency of the trigger selection as 
the software trigger changed during the data taking in 2012. The early 2012 setup is not 
modelled in the simulated samples, and so it is assumed that the efficiency of this stage of 
the trigger selection is the same for both periods. The uncertainty due to this assumption 
is taken to be 5%, which is consistent with differences in the efficiencies measured in 
analyses of similar decays. Since this data comprises ~ 20% of the total data sample, the 
overall uncertainty on the trigger efficiency is 1%. 

The reconstruction efficiency of the photons cannot be determined using the simulated 
samples, and so a data-driven method is used to correct for this efficiency 42 . The 
correction factor, C'.y, is found in bins of photon pt by comparing the relative yields of 
reconstructed B+ — )■ J/r|;K*+(— )■ K+7t°) and B+ — )■ J/r|;K+ decays. Since the A°— )■ Arj' 
(rj' —>■ Tt+Tfy) decay contains the same number of photons as the normalisation channel, this 
correction cancels in the ratio of efficiencies. However, for the A°—)■ Arj' (ri^— 7f^7r”r|) and 
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A°—)■ Arj (t]—)■ decays, there is an extra photon in the signal channel compared 

to the normalisation channel. The correction factor is therefore applied to these channels, 
and an uncertainty is introduced due to the limited size of the data sample used to measure 

it. 

The systematic uncertainty of the PID selection is estimated by the difference between 
the efficiency estimates obtained in data from calibration samples, and the efficiency 
extracted from the simulation after reweighting the PID variables to match the respective 
distribution observed in data. 

Finally, there is an uncertainty due to the procedure for handling events with more 
than one candidate. The efficiency for selecting the best candidate using the procedure 
described earlier, and the fraction of events which contain more that one candidate are 
both estimated using simulated samples of signal decays. The estimated uncertainty on 
each of these quantities is around 10%, which is chosen as a conservative uncertainty. This 
introduces an overall systematic uncertainty of 1.4% to 2.2% depending on the signal 
channel. 


6 Confidence Intervals 


The number of signal events observed in each of the signal decays can be used to place a 
limit on the ratio of branching fractions of the signal decay with respect to the normalisation 


channel. The unihed approach method presented in Ref. 43 is used to place a limit on 


the branching fraction. This method constructs confidence intervals based on a likelihood 
ratio method using the probability of observing N^hs signal events for a given branching 
fraction. For this analysis the probability is assumed to follow a Gaussian distribution 
with a resolution of a = where cTsyst is the systematic uncertainty described 

above and cxstat is the statistical uncertainty on the number of events observed. 

The 68% and 90% conhdence level (CL) intervals are obtained for the A°—;■ Ap and 
A°—)■ Ap' decays respectively, by combining the likelihoods of each category. The weighted 
average of the observed ratio of branching fractions is calculated from the number of events 
observed in each signal channel with their uncorrelated systematic uncertainties and the 
values of a in Table and this is used to construct the conhdence intervals. The combined 
L and D categories are used to construct conhdence intervals for the Aj( —)■ Ap decay. 
The likelihoods from the A°—)■ Ap' (p'—?■ Tr+Ti'y) and A{^—Ap' (p'—)■ 7t+7t“p) decays are 
then combined, with the L and D combined, to give a limit on the branching fraction 
for the A^—)■ Ap' decay. The limit on the ratio of branching fraction with respect to the 
normalisation channel is found to be 


g(Ag^Ap') 

KV) 


< 0.047 at 90% CL, 


for the A°—)■ Ap' decay, and for the A°—;■ Ap decay the 68% CL intervals are 

i3(B0^KV) 
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Multiplying by the known value of B (B° —)■ K^rj') gives a limit on the branching fraction 
of the A° —)■ Ari' decay 

i3(A0^ Ati') < 3.1 X 10"® at 90% CL, 
and the 68% CL intervals for the A° —)■ Ar| decay are 

i3(A0^ATi) = (9.31^1) xlO-^ 


7 Conclusions 

A search is performed for the A° —)• Ap' and A° —)• Ap decays in the full dataset recorded 
by the LHCb experiment during 2011 and 2012, corresponding to an integrated luminosity 
of 3fb~^. No signihcant signal is observed above background for the A°—)■ Ap' decay, and 
some evidence is seen for the A{^ —)■ Ap at the level of 3 ct. The B° —)■ K°p' decay is used 
as a normalisation channel, so that a limit is placed on the ratio of A°—)• Ap^'^ branching 
fractions with respect to the B° —)■ K'^p' branching fraction using the unihed approach. 
With the known value of the B*^ —)■ K°p' branching fraction, the upper limit on the absolute 
branching fraction of the A°—)■ Ap' decay is i3(A°—)■ Ap') < 3.1 x 10“® at 90% CL. The 
branching fraction of the A° —)■ Ap decay is S(A0^Ap) = (9.3;^J)xl0-6. 

These values can be compared with the branching fractions calculated in Ref. (^ , 
and given in Sec. The predicted branching fractions depend strongly on the method 
used to calculate the hadronic form factors, and on the parameters used in the calculation, 
as discussed earlier. Our results favour the branching fractions calculated using the pole 
model to estimate the hadronic form factors. In addition, our results are inconsistent with 
the prediction for B (A° —)• Ap') obtained by neglecting the anomalous contribution to the 
decay amplitude, indicating that a gluonic component of the p' wavefunction should be 
present. 
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